Reentrant topological phase transitions in a disordered spinless superconducting wire 
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In a one-dimensional spinless p-wsLwe superconductor with coherence length ^, disorder induces a phase 
transition between a topologically nontrivial phase and a trivial insulating phase at the critical mean free path 
/ = ^/2. Here, we show that a multichannel spinless p-wave superconductor goes through an alternation of 
topologically trivial and nontrivial phases upon increasing the disorder strength, the number of phase transitions 
being equal to the channel number A'. The last phase transition, from a nontrivial phase into the trivial phase, 
takes place at a mean free path / = ^/(N + 1), parametrically smaller than the critical mean free path in one 
dimension. Our result is valid in the limit that the wire width W is much smaller than the superconducting 
coherence length ^. 
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In one dimension, spinless superconductors appear in two 
topologically distinct phases. In one of these phases, usually 
referred to as the "trivial phase" the excitation spectrum is 
adiabatically connected to the ionic insulator. The other phase 
is "topologically nontrivial". Topologically protected zero- 
energy bound states appear at junctions between the trivial 
and nontrivial phases |[Tl|2l. These bound states are particle- 
hole symmetric and two of these combine to form a single 
fermionic excitation, which is why they are referred to as 
"Majorana bound states" |T, 4|. Interest in these systems has 
peaked after recent proposals to construct topological super- 
conductors out of hybrid structures involving standard BCS 
superconductors and semiconductors ||5]|6l and reports of their 
subsequent experimental realization QUI- 

The Pauli principle enforces that spinless superconducting 
correlations are odd in momentum. In a one-dimensional set- 
ting, this means that they must be of /j-wave type. Unlike for 
s-wave superconductors, where the Anderson theorem pro- 
tects the superconducting correlations against impurity scat- 
tering [9J, backscattering by impurities suppresses p-wave su- 
perconducting order. In a one-dimensional wire any small 
amount of disorder already leads to subgap states at arbitrarily 
low energies, but it takes a finite amount of disorder to drive 
the system from the nontrivial superconducting phase into the 
trivial phase ifTOl [TTl . For short-range disorder with normal- 
state mean free path the transition between these phases 
takes place if IfTTl 

(1) 

where ^ is the superconductor coherence length. Here and 
below we assume that the superconductivity is weak, ^ much 
larger than Fermi wavelength Ap. 

The one-dimensional description applies only if the system 
width W does not exceed the Fermi wavelength Ap. If W > Ap, 
the normal-state has - int {2W/Ap) > 1 propagating chan- 
nels at the Fermi level, and without disorder the topologi- 



cally nontrivial phase exists if is odd, but not if is even 
lfT2l - nTl . Numerical simulations and weakly-disordered per- 
turbation theory indicate that the topological phases are stable 
against weak disorder in the multichannel case, too 11811191 . 

It is the purpose of this letter to provide an analytical the- 
ory of the effect of disorder on the topological phase in the 
A^-channel p-wave superconductor. Our main result is that 
increasing the disorder strength drives the system through a 
sequence of A^ topological phase transitions, taking place at 



In particular, a topologically nontrivial phase persists for dis- 
order strengths up to / = ^/(N + 1), significantly larger than 
the critical disorder strength ([T]l at which the topological phase 
transition takes place in one dimension. Our analytical theory, 
as well as the precise location of the phase transitions given 
in Eq. Q, is valid in the limit of thin wires, width W <*c ^. 
We have verified numerically that the alternation of topologi- 
cal phases persists for wire widths up to W ~ ^. Note that the 
existence of A^ phase transitions is consistent with the known 
results for the weak disorder limit I — > oo (nontrivial phase if 
A^ is odd, and trivial phase if A^ is even), as well as the strong 
disorder limit / i (system is in the trivial phase). 

For the derivation of Eq. (j2|i we consider a spinless p-wave 
superconducting wire of length L, width W, and chemical po- 
tential p > coupled to ideal normal-metal leads at its two 
ends. The Bogoliubov-de Gennes Hamiltonian of the system 
reads 

"^{h^ ^(-^.y) - juj cr, + ^ {A;, p,} cr_, + A'^pyCr,., (3) 

where ct j, ct,,, and cr, are Pauli matrices in particle-hole space. 
The Hamiltonian ([3]) has particle-hole symmetry, o-xHa-^ - 
-H*, which places it in the (Altland-Zirnbauer) symmetry 
class D |20|. The superconductor occupies the volume < 



2 



X < L, and the two leads are at x < and x > L respectively, 
see Fig. [T| (inset). The superconducting order parameters 
and A,, are nonzero for < x < L only. The superconducting 
coherence length is 



(4) 



Although A\ - A', for an isotropic superconductor, we have 
chosen to use different symbols in order to underline the very 
different roles of these two parameters in the calculation that 
follows. We assume that the superconductivity is proximity- 
induced, so that we can treat A\ and A', as externally-imposed 
parameters without self-consistency requirements. The im- 
purity potential has zero average and short-range fluctuations 
described by the Gaussian white noise correlator 



{V{x,ymx' ,y')) = y6{x - x')5{y - y') 



(5) 



and is zero in the leads. 

We determine the topological phase from the zero-energy 
reflection matrix r of the superconducting wire ||2T]| . In the 
particle-hole notation, the reflection matrix r for quasiparti- 
cles incident from the left takes the form 



With the chiral symmetry present it is possible to ex- 
press the zero-energy reflection matrix r in terms of the sys- 
tem's normal-state scattering matrix at a slightly renormal- 
ized chemical potential jl |23|. To this end, we first rotate 
the Hamiltonian Q to the Majorana basis 

(9) 



At zero energy, the eigenvalue equation for H consists of two 
decoupled equations describing particles that are exposed to 
an imaginary "gauge field" of magnitude ^/^ and pointing 
in opposite directions for the two equations ll24l l25l . This 
"gauge field" may be transformed away by the (non-unitary) 
transformation 



tl/(x,y) ^±{x,y) = < 



ip{x,y) 

4i{x,y)e~ 

\p{x,y)e^^l^ 



if X < 0, 
if < X < L, 
if X > L. 



(10) 



The wavefunctions ij/-t(x,y) satisfy the standard Schrodinger 
equation for the zero-energy wavefunction of a disordered 
wire. 



^ee ^eh 
'"he '"hh 



(6) 



and 



where particle-hole symmetry imposes that rhh = 
The = r*j^ at zero energy. Following Fulga et ai, the topolog- 
ical phase can be calculated from the determinant Q = det r 
ETl : The topologically nontrivial phase has 2 = -1, whereas 
the topologically trivial phase has 2=1. (Note that particle- 
hole symmetry requires det r to be real; As no extended quasi- 
particle states exist in the superconductor away from the crit- 
ical points, r must be unitary and hence | det r\ - 1.) 

In the thin-wire limit <k ^ the transverse pairing AJ, may 
be treated perturbatively 1 17|. Without the transverse pairing, 
the Bogoliubov-de Gennes Hamiltonian H has an additional 
chiral symmetry a-yHcTy - -H |22|, which places it in the 
symmetry class BD I. With the chiral symmetry, the topolog- 
ical superconducting phases are characterized by an integer 
number Qchirai- The topological quantum number Q is related 

to gchiral as 



G = (-i) 



Qchiml 



(7) 



The absolute value IQchirail can be interpreted as the number 
of Majorana bound states at the end of the wire, when the 
normal metal leads are replaced by insulating ends ETl (see 
also App. 2). The quantum number Qchirai can be calculated 
from the zero-energy reflection matrix r as i2H 



chiral 



-i lim tiTeh- 

L— »oo 



(8) 



The limit L — > oo is taken in order to ensure that the reflection 
matrix r is unitary 121]. The chiral symmetry implies that r^h 
is an antihermitian matrix, reh = ~'"eh' so that Qchirai is real. 



2m 



-fi + V\tfr^{x,y)^0, 



(11) 



where ft = fi + fiA'Jl^. Transforming back to the basis of the 
original Hamiltonian ^ allows us to express the reflection 
matrix r in terms of the reflection matrix f (for particles inci- 
dent from the left) and the transmission matrix F (for particles 
incident from the right) of the normal-state scattering problem 



specified by Eq. ( 1 1 



[1 +t'rsinh^(L/^)] 



reh = /sinh(L/^)cosh(L/f)[l + t'r sinh^(L/^)]"'f'r"f. (12) 

Returning to Eq. ([8]l, we find that topological number Qchiiai 
can be expressed as a sum over the eigenvalues f„(L) of f'P' , 



chiral 



lim y 



n=l 



sinh(L/^) cosh(L/^)f„(L) 
l+f„(L)smh\L/^) 



(13) 



The asymptotic probability distribution of the eigenvalues 
f„ in the limit of large L is well studied in the literature f26'|. 
The result is best parameterized in terms of the "Lyapunov 
exponents" f„ - cosh"^(x„L), which are self-averaging for 
large L, with mean 



{x„} 



{N+l)l 



1,2, 



(14) 



and small fluctuations of order 1 / y/iN + 1)IL. The mean free 
path is / - Afh^viai^ l{2y), with a numerical factor de- 



pending on the width W. Substituting this result into Eq. (13 1, 
we find that 



chiral 



N 
n=l 



(N + 1)1 



(15) 
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where 0(z) = if z < and 1 otherwise. Hence, upon in- 
creasing the disorder strength, the topological quantum num- 
ber 2chirai Stepwise decreases from Qchirai = in the limit of 
zero disorder to Qchhai = in the strong disorder limit. The 
transitions take place at the critical disorder strengths of Eq. 
Q. The topological quantum number Q is given by Eq. (|7|. 

When the transverse coupling proportional to A', is taken 
into account, the chiral symmetry is broken, and the topo- 
logical quantum number QcWiai is no longer meaningful. The 
quantum number Q remains well defined, however. Since the 
effect of A,, is small if W <s; ^ 1 17 1, the value of Q remains un- 
changed upon inclusion of the transverse coupling. Upon in- 
creasing the disorder strength, we therefore expect alternation 
between topolgical trivial {Q - 1) and nontrivial {Q - -1) 
phases until, in the limit of strong disorder, / < ^/(N + 1), the 
system remains in the trivial state. As long as W « ^, the 
transition points should exhibit only a weak dependence on 
the transverse coupling A[.. 

Alternatively (and equivalently), for a superconductor wire 
with hard-wall ends, the transverse coupling pairwise gaps out 
the 2chiiai Majorana bound states at each end of the supercon- 
ducting wire, leaving behind a single Majorana state if and 
only if Qchirai is odd. Since Qchiiai decreases stepwise from 
to zero upon increasing the disorder strength, the number of 
Majorana bound states at the end of the wire with A[, taken 
into account alternates between and 1, the transitions taking 
place precisely at the disorder strengths given in Eq. Q. Since 
the presence or absence of a single Majorana fermion is topo- 
logically protected, inclusion of the transverse coupling A[,py 
for sufficiently small W/^ does not affect these transitions or 
the transition points. 

For broader wires, W ~ ^, the transverse coupling can not 
be treated perturb atively, and the results for the chiral limit 
A', — > at best have qualitative validity if the transverse cou- 
pling is included. In this respect, we note that the topological 
phase transitions no longer take place at weak disorder I » Ap 
if W ~ ^. This can be seen from Eq. (|2| upon substituting 
W ~ (N + l)Ap, which gives / ~ nAp. Equation ( 14 1, which 
was essential for establishing the transition points, is derived 
under the assumption of weak disorder, / » Ap 82611 . and no 
longer has quantitative validity if this condition is violated. 

In order to further support our conclusions and to investi- 
gate the regime W ~ ^,we have performed numerical simula- 
tions of the Bogoliubov-de Gennes Hamiltonian We cal- 
culate the reflection matrix r by concatenating short segments 
of length 6L <s: Ap. We refer to Ref. 27 for a description of the 
numerical method. The scattering matrix of a short segment 
is calculated to lowest order Born approximation. For techni- 
cal reasons the numerical data were obtained by varying the 
magnitude of the superconducting parameters A'^ and A,, and 
keeping a fixed mean free path /. 

First, we verify our analytical results in the chiral limit, 
A', - 0. Figure |Tj shows the topological number Qchiiai as a 
function of the ratio ^/l for a wire with N - 9 channels. The 
figure clearly shows the stepwise decrease of Qchiiai upon in- 
creasing the disorder strength in comparison to the supercon- 




^/(N-t-l)l 

FIG. 1. (Color online.) Topological number gcWrai for a wire with 
width IWjAf = 9.5 such that the channel number is A' = 9, as 
a function of the ratio ^// of disorder strength and induced super- 
conductivity. Data shown are for a single disorder realization with 
Ap/l = 0.011 and wire length L/l ~ 2100. The red curve shows 
the analytical prediction ijTSj and the blue one the numerical data. 
Inset: Schematic picture of a disordered superconducting wire with 
two ideal normal-metal leads. 



ducting order For coherence lengths of the order of the local- 
ization length {N + 1)1 the transition points closely follow the 
theoretical prediction (|2]i. We attribute the quantitative devi- 
ation of the transitions at large Qchiiah when ^ ~ I, where the 
relevant Lyapunov exponent x„ is comparable to the inverse 



mean free path, to a failure of the estimate ( 14 1 in this regime 
I26j. 

The effect of the transverse coupling on the phase transi- 
tions is shown in Fig. |2] Both panels of Fig. |2]show the value 
of Q - detr as a function of (A^ -i- 1)//^ and of A',,/A^ at a 
fixed realization of the random potential V (i.e., at a fixed 
value of the mean free path /). The top panel of this figure 
shows representative numerial data for a weakly disordered 
wire (Ap/l = O.OIl), where all transitions take place within 
the limit W <^ ^. As expected, the sequence of topologi- 
cally trivial and nontrivial phases does not significantly de- 
pend on the transverse coupling A[. in this case. The bottom 
panel shows data for strong disorder (Ap/l = 0.43), where the 
condition W <s; ^ is no longer satisfied for small values of 
^/l. The disorder strength is chosen such that the condition 
W - ^ is met roughly at the 6th transition. For the necessarily 
finite wire lengths L in the numerical simulations, finite-size 
effects lead to a blurring of the topological phase transitions. 
The occurence of values of det r different from - 1 or 1 signals 
a breakdown of the insulating behavior of the superconduc- 
tor This behavior is consistent with Ref. 1281 . where it was 
found that for large a spinless superconducting wire enters 
a quasi-critical region with algebraic instead of exponentially 
decaying transmission |29|, which persists up to wire lengths 
L much larger than the normal-state localization length and 
out of range of our numerical simulations. It is also consis- 
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FIG. 2. (Color online.) Phase diagrams showing the topological 
number Q = det(r) of a wire with 2W/Af = 9.5 as a function of the 
ratio //^ and the transverse coupling A', in a spinless superconducting 
wire witih N = 9 channels in the limit of weak and strong disorder 
(panel a) and b), resp.). The disorder strength is Ap/l = 0.0 11 (panel 
a)) and A^/l = 0.43 (panel b)). The data were obtained for a wire 
length L/l = 150 (panel a)) and L/l ~ 770 (panel b)). 



tent with the approach of the two-dimensional limit, in which 
the one-dimensional thermal insulator transitions into a two- 
dimensional thermal metal ll30l . 

In conclusion, we investigated the effect of disorder on the 
topological phase in a multichannel p-wave superconducting 
wire. From an analytical study in the limit of thin wires 
W <s; ^, we derived a series of topological phases transitions 
in which the system alternates between trivial and nontrivial 
phases. A numerical analysis shows that this holds true also 
for thicker wire W < ^ 

We gratefully acknowledge discussions with Michael Wim- 
mer and Falko Pientka. This work is supported by the Alexan- 
der von Humboldt Foundation and by the funds of the Erdal 
Inonu chair 



Appendix 1: Full scattering matrix. 



used the mapping to the normal scattering problem ( 1 1 



In the main text, we 
to 



compute the left reflection matrix of the superconductor. The 
same mapping can be used to compute the full scattering ma- 
trix. We here summarize the main resuls. 

The scattering matrix of the normal metal wire may be writ- 



^4; 



U 
V 



- Vi -f 



yT 



(16) 



where U and V are unitary matrices and f = 
diag (fi, f2, . . . , f„). For the scattering matrix S of the 
superconductor we then find 



r t' 
t r' 



(17) 





U* 





0^ 





V* 



PPT 

-p -p -f 

T f p 
-f T P 



f ^ 

T 
-P 
P I 











t/-'- 









yT 





0^ 







where we abbreviated 



[l+fzl 



p^i\l+ fz^_ 



r'vr 
1" 



Z+Z-T 



[l+fz^J 



t = ;■ [l -I- fz^] ' z- V^O"^-^, 



(18) 



with z+ - cosh(L/^) and Z- - sinh(L/^). 

Appendix 2: Number of Majorana fermions. To connect 
the topological quantum number Qchirai to the number of Ma- 
jorana fermions for a superconducting wire with hard-wall 
ends, we recall that the number of Majorana fermions is equal 
to the dimension of the nullspace of 1 + r. Making use of 



our explicit solution (17 1 for the reflection matrix r in terms 
of diagonal matrices p and p, we see that the elements of the 
diagonal matrix p are or 1 in the limit L — » oo, and that 
p - i(\ - p). We may then factorize 1 + r as 



\+r^U 



UW*-p i{l-p) 
-i(l-p) U'^U-p 



U 



-iU^ 



ip 



(19) 



ipi\p 

From this expression one immediately concludes that the di- 
mension of the nullspace of 1 -i- r is the number of unit eigen- 
values of -ip, which is equal to Qchiiai by Eq. ([8|. 

Appendix 3: Mean free path The mean free path in the 
normal state is calculated from the reflection matrix for a short 
segment of length dL as tr^ree?"fie) = NdL/l, which is calcu- 
lated in the first order Born approximation. The expression 
for I given in the main text contains a factor that explicitly 
reads 



On - N 



3 -1 
^Y^i^Wkp/nf-n^y 



n=l 



+2 J] yjaWkF/nf - n^rHiWkF/nf - n^Y 



(20) 
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